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o 

s ! ' Abstract. We carry out some computations of vector valued Siegel 

(-H , modular forms of degree two, weight (k, 2) and level one. Our approach 

is based on Satoh's description of the module of vector-valued Siegel 
modular forms of weight (k, 2) and an explicit description of the Hecke 

i~^ . action on Fourier expansions. We highlight three experimental results: 

(1) we identify a rational eigenform in a three dimensional space of cusp 
forms, (2) we observe that non-cuspidal eigenforms of level one are not 
always rational and (3) we verify a number of cases of conjectures about 

S congruences between classical modular forms and Siegel modular forms. 

-i— > 
S ■ 

1. Introduction 

Computations of modular forms in general and Siegel modular forms in 
^ \ particular are of great current interest. Recent computations of Siegel mod- 



in 



ular forms on the paramodular group by Poor and Yuen 20J] have led to the 
careful formulation by Brumer and Kramer [3J of the Paramodular Conjec- 
ture, a natural generalization of Taniyama-Shimura. Historically, computa- 
tions of scalar-valued Siegel modular forms in the 1970s by Kurokawa [161 ] 
f-^ \ led to the discovery of the Saito-Kurokawa lift, a construction whose gen 



£N) ■ eralizations are still studied. In the 1990s, computations by Skoruppa 28| 

revealed some striking properties that some Siegel modular forms possess 
(namely, there are rational eigenforms of weights 24 and 26 in level 1 that 
span a two-dimensional space of cusp forms). These properties have yet to 
j_j \ be explained. This paper is in the same spirit. 

We carry out the first systematic computations of spaces of vector- valued 
Sieg el modular forms of degree two and of weight (k, 2). We do this in Sage 
[25[ using a package co-authored by the second author, Raum, Skoruppa 
and Tornaria [22]. We observe some new phenomena (see Propositions 13.21 
and 13 . 3 j) and check that the eigenforms we compute satisfy the Ramanujan- 
Petersson bound (see Proposition 13. ip . 

We also verify two interesting conjectures on congruences: the first, due to 
Harder, has been previously verified by Faber and van der Geer [S|; the other, 
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due to Bergstrom, Faber, van der Geer and Harder, has been previously 
verified by Dummigan [21] . Our approach to verifying these conjectures is to 
compute Hecke eigenvalues of Siegel modular forms in as direct a manner as 
possible, using Satoh's concrete description of Siegel modular forms of weight 
(k, 2). This is a very different approach than the one taken by Faber and 
van der Geer and we verify cases that they do not (and vice versa). After 
determining a basis of eigenforms for the space, we use explicit formulas 
for the Hecke action on Fourier expansions to extract the Hecke eigenvalues. 
Our main results in this direction are Theorem l4.3l and l4.7l which summarize 
the cases of the two conjectures that we have verified. 

We do not describe the implementation of Siegel modular forms used to 



carry out these computations but we refer the interested reader to 2l| for 
such a description. The point of this paper is that we were actually able to 
carry out such computations, have made some new observations based on 
these computations and have made our data publicly available [24|. 

2. Vector- valued Siegel modular forms of weight (k, 2) 

We recall the definition of a Siegel modular form of degree two. We 
consider the full Siegel modular group T^ 2 ' given by 

r^ := Sp(4, Z) = JM G M(4, Z) : t M (-i 2 h ) M =(-i 2 l2 )}- 

Let 

M (2) := {Z € M(2,C) : l Z = Z,Im(Z) > 0} 

be the Siegel upper half space of degree 2. For a nonnegative integer j, the 
space C[X, Y]j of homogeneous polynomials of degree j has a GL2-action 
given by 

(A,p)^A-p:=p((X,Y)A). (2.1) 

Definition 2.1. Let k,j be nonnegative integers. A Siegel modular form of 
degree 2 and weight (k, j) is a complex analytic function F : H' 2 ) — > C[X, Yh 
such that 

F(gZ) := F((AZ + B){CZ + D)' 1 ) = det(CZ + D) k (CZ + D) ■ F(Z) 
for all g= ($g) G r( 2 ). 



Remark 2.2. The action in (12. lj) is a concrete realization of the symmetric 
power representation Sym^ of GL2 . Definition 12.11 is a concrete description 
of Siegel modular forms with values in the representation space det (8> Sym 5 
of GL2 . We made these choices in our implementation because it made the 
multiplication of vector- valued Siegel modular forms easier to implement. 

(2) 
The space of all such functions is denoted M, ■ , where we suppress j if it 

is 0. If j is positive F is called vector-valued, otherwise it is called scalar- 

(2") (2^ 

valued. We write M* := ® fe M k for the ring of (scalar- valued) Siegel 
modular forms of degree 2. 
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Let 

Q-= [f = [ a , b, c] : a,b,c£ Z, b 2 - Aac < 0, a > 0} 

where [a, b, c] corresponds to the quadratic form aX 2 + bXY + cY 2 . 
A Siegel modular form F has a Fourier expansion of the form 

F(Z)= Y, C F (f)e(aT + bz + CT'). 

f=[a,b,c]£Q 

Here Z := {J Z Z T .) (t,t' G H* 1 ) and z G C), e(x) = e 2nix , the trace of a 
matrix A is denoted by tr A. The form F is called a cusp /orm if its Fourier 
expansion is supported on positive-definite elements of Q. The subspace of 
cusp forms is denoted Si ■ . 

(2) 

The ring of all vector- valued Siegel modular forms (J)fc 7 - M^ is not finitely 
generated. For this reason the symmetric power j is usually fixed. The 
resulting module is finitely generated over M* . We focus exclusively on 
weight (k,2), where we have a very concrete description of these spaces 
thanks to work of Satoh. 



2.1. Satoh's Theorem. The Satoh bracket is a special case of the general 



(2) (2) 

(k, 2). Suppose F E M^ and G G My are two scalar- valued Siegel modular 



Rankin-Cohen bracket construction. Satoh 26] examined the case of weight 

(Jfc, 2). Suppose F G M^ ] and G G M$ 
forms. We define the Satoh bracket by 

[F,G]2 = ^~ {^Gd z F - yFd z G^j G M^, 2 , 
where d z = (J^Jf^)- 

(2) 

In the same paper, Satoh showed that fe M^ J, is generated by elements 
all of which can be expressed in terms of Satoh brackets. More precisely, he 
showed that 

M® =[E 4 , E 6 ] 2 ■ Mi 2 } 10 [E 4 , Xl0 ] 2 ■ Mf_ } 14 e 

[E-4, Xnh ■ MJ® 16 [E G , X ioh ■ C[E 6 , xio, Xulfc-w© 

[E(5,X12]2 ■ C[E 6 ,XlO,Xl2]k-18 © [XlO, Xl2]2 ■ C[X10, Xl2]k-22- 

Here the forms E^, Eq,xwi X12 are the generators of the ring of scalar- valued 
Siegel modular forms described by Igusa [12]]. By CLAi, . . . ,A n ]k we mean 
the module of weight k modular forms that can be expressed in terms of 
generators A\, . . . , A n . 

A basis for the space Mf, 2 was computed via a Sage [25(] implementation 



in [22| of an algorithm found in |21l |. 
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2.2. Hecke Operators. As our interest is in computing Hecke eigenforms, 
we need to describe how one computes the Hecke action. We give formulas 
for the image of a Siegel modular form of weight (k, 2) under the operator 
T(p ). The Hecke operators are multiplicative and so it suffices to under- 



stand the image for these operators. The formulas can be found in 11] but 
we present them here for completeness. 

Let F be a Siegel modular form as above and let the image of F under 
T{p ) have coefficients C'([a,b,c\). Then 



C'([a,b,c})= Yl p^+^-^x 

au bu cu 



au^OipP+l) 
b u =c u =0(p~ l ) 



pP+1 ' pi ' p-1-P 



(2.2) 



where 



R{y) is a complete set of representatives for SL(2, Z)/iq {jpP) where 

Tg (jr) is the congruence subgroup of SL(2, Z) of level p 13 ; 
for f = [a,b,c], [a u ,b u ,c u ] = f u :=f((X,Y) t U)- 

l 

p^ 

the • is given by the action defined in (|2.1|) 



<k,fi = I „/3 



We denote the Hecke eigenvalue of a Siegel modular form F under the opera- 
tor T(p s ) by X p s(F). If the space S k 2 has dimension d, the Hecke eigenvalues 
of F are algebraic numbers of degree at most d. The field that contains the 
Hecke eigenvalues of F is denoted Qf . 

2.3. Computing Hecke eigenforms. Fix a space of Siegel modular forms 
of weight (k, 2) with basis {F\, . . . , F n }, obtained as algebraic combinations 
of the Igusa generators and Satoh brackets. Because the Hecke operators 
are a commuting family of linear operators, there is a basis {G\, . . . , G n } for 
the space consisting entirely of simultaneous eigenforms. 

The forms G{ are determined computationally as follows. First, deter- 
mine the matrix representation for the Hecke operator T(2) by computing 
the image under T(2) of each basis element F^. Build a matrix TV that is 
invertible and whose jth row consists of coefficients of Fj at certain indices 
Qi, ■ ■ ■ i Qn- To ensure that N is invertible we pick the indices one at a time, 
making sure that each choice of index Qi increases the rank of N. We then 
construct a matrix M whose j'th row consists of coefficients of the image of 
Fj under T(2) indexed by Q\, . . . , Q n . Then the matrix representation of 
r(2) is MN- 1 . 

We compute the Hecke eigenforms using T(2) and we express them in 
terms of the basis {F±, . . . ,F n }. We compute the Hecke eigenvalues A p <s 
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by computing these expressions to high precision and then computing their 
image under the Hecke operator T(p ) as in (j2.2p . 

2.4. Hecke eigenvalues, Satake parameters and symmetric poly- 
nomials. Fix a prime p. The Satake isomorphism f2 is a map between 
the local-at-p Hecke algebra 7~L P associated to V and a polynomial ring 
Q[xq,xi,x 2 ] W2 invariant under the action of the Weyl group. A matrix 
representation of Q, can be found in 14j, |23J. We summarize the relevant 
results here. 



Consider 



T{jp) 



Up 2 



V 

/I 

V 



p 



p 



\ 



P) 



v) 



Up 2 ) 



T2{P 2 ) 



V 

(p 

{ 



\ 



p 2 ) 



pj 



The images under fi of these operators are: 

Q(T(p)) = X XiX 2 + XqX\ + X X 2 + x 

tt(T (p 2 )) 



n(T x ( P 2 )) 

n(T 2 (p 2 ) 



2p-2 
2 -l 



>2 + ^</>l + 



/' 



^92 



>2 + - p (f>l 



where 



222, 22, 22, 2 
XqX]^X2 i~ XqX^ -}- XqX2 i~ Xq 

XqX^X 1 } \ XqX\X2 ~t~ XqX\ -\- XqX2 



4>o = 

01 = 

02 = X5X1X2. 

Fix a Siegel Hecke eigenform F. It can be shown [l|, p. 165] that for any 
p there exists a triple (a^ , af^ , a^ p ) € (C x ) 3 /W / 2 with the property that 
Q(T)\ f = Xt(F), the eigenvalue of F with respect to the Hecke oper- 

ator T. The numbers a are called the Satake parameters of -F at p. 

We will make use of the following way of expressing T(p) 2 in terms of the 
operators T,(p 2 ): 

Theorem 2.3 (@|). T{p) 2 = T {p 2 ) + (p + ljT^p 2 ) + {p 2 + l)(p + l)T 2 {p 2 ). 

3. Computational and Experimental Results 

We carry out the computations of particular eigenforms in the following 
way. Satoh's theorem as described above in Section 12.11 gives a recipe for 
computing a basis of vector- valued Siegel modular forms of weight (k,2). 
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In particular, following 28J] we can compute the Fourier expansions of the 
Igusa generators E^Eq, X10jX12- This is done via an explicit map from 
elliptic modular forms to Siegel modular forms. For each generator we easily 
computed the part of its Fourier expansion which is supported on positive 
definite quadratic forms up to discriminant 3000 and on singular quadratic 
forms [0, 0, c] where < c < 750. 

Using these four Igusa generators we determine a basis for the space of 
weight (k, 2) as prescribed by Satoh's theorem: we compute a basis for the 
modules M fc _i , M fc _i 4 , M k _ 16 and C[£ 6 ,xio, X12U-I6, C[£ 6 ,xio,Xi2]fc-i8, 
< C[xiO;Xl2]fc-22- We then form a basis for the space of weight (k,2) by 
multiplying each basis above by the Satoh bracket that corresponds to it 
in Satoh's theorem and end up with eigenforms following the procedure 
described in Section \2. 31 One might pause at the idea of multiplying vector- 
valued Siegel modular forms but this is precisely the reason why we defined 
the coefficients of vector-valued Siegel modular forms to be homogeneous 
polynomials (see Remark I2.2p . 

For example, we find that a basis for the space of weight (16, 2) is given 
by 



n — VP \VP Pl 173820100608 r p .. 1 , 1800409600 r F ,, 1 
<-n - m [&4,&6 1 2 1557539 — [^A,XV2\2 + 1557539 L^6 , X10J2 

G 2 = [E^Xuh + (s064 Q + "uM-^jXlofe 

Gi = [i?4,Xl2]2 + (§064" + l2")[-^6,Xlo]2 



where a is a root of x 2 + 58752x + 858931200 and a is its conjugate. The 

form G\ is non-cuspidal (probably Eisenstein) and the other two forms in 

the basis are cuspidal. 

(2) 
The Hecke eigenvalues that appear in the space SI 2 tend to have the 

largest possible degree, namely the dimension d of the space. For example 
the forms G2 and G2 above have Hecke eigenvalues in a quadratic field. 
There are however (very surprisingly!) counterexamples to this; this is anal- 
ogous to what happens in the scalar- valued spaces 5^ and S^q , see [281 ] . 
but in stark contrast to the situation in degree one (as predicted by Maeda's 
conjecture). 

Consider the weight (20, 2). Let K = Q[a] be the quadratic number field 

with minimal polynomial x 2 — 780288x + 121332695040, and let a denote 

(2) 
the conjugate of a in K. The space 5^2 is three-dimensional, and a basis 
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of Hecke eigenforms is given by: 
Hi = Xio[E4,E 6 ] 2 — j^Eq[E4,xio] 

H 2 = Xw[E4,Eq} 2 + ( 12241152 ° ~~ 15939 ) E e[ E i, Xlo]2 + 



(354520° + m) E 4 E ^ xi2J2 

H3 = Xlo[-E , 4,-E , 6]2 + ( i2241152 ° ~ T5939)- E '6[-E , 4,Xlo]2 + 

(363520° + m) E *[ E *> ^ 12 ]2- 

Note that the first has rational eigenvalues, and the second and third 
have conjugate quadratic eigenvalues. We checked that each of these forms 
satisfies the Ramanujan-Petersson conjecture. 

Proposition 3.1. For all k satisfying 14 < k < 30, the Hecke eigenforms 

(2) 
in Su2 satisfy the Ramanujan-Petersson conjecture at p = 2,3,5. More 

(2) 
precisely, let F 6 Sj,^ be a Hecke eigenform with eigenvalues \ p , X p 2 and 

consider the polynomial 

X 4 - A P X 3 + (Aj - X p2 - P 2k - 2 ) X 2 - v 2k - x \x + p 4k ~ 2 . 

Then all roots z € C of this polynomial satisfy \z\ = p( 2fc +J -3 )/ 2 . 

We also looked at the level 1 elliptic modular forms with Hecke eigenvalues 
in quadratic fields and these fields are different than K. This indicates that 

H2 and H% are unlikely to be lifts. Therefore the naive generalization of 

(2) 
Maeda's conjecture does not hold for Sio- We remark that in all other 

weights for which we have carried out computations, the naive generalization 

does indeed hold: 

Proposition 3.2. Let k G {14,16,18,22,24,26,28,30}. Then the charac- 
teristic polynomial of the Hecke operator T(2) acting on S^.2 is irreducible 
over Q. If k = 20, the characteristic polynomial of the Hecke operator T(2) 
decomposes over Q into a linear factor and a quadratic factor. 

We also make note of another interesting computational phenomenon that 
merits further investigation. It is interesting to us because an analogous 
phenomenon does not happen in the scalar-valued case. In the scalar-valued 
case of level 1, there are four kinds of modular forms of even weight: Eisen- 
stein series, Klingen-Eisenstein series, Saito-Kurokawa lifts and cusp forms 
that are not lifts. The first two are not cuspidal and always have rational 
coefficients. Compare this fact to the following proposition: 

Proposition 3.3. Let k € {22,26,28,30}. The space of modular forms of 
weight (k,2) that are not cusp forms is two dimensional but consists of a 
single Galois orbit. 

The data that make up the proofs of these propositions can be found at 



24] 
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4. Verification of some conjectural congruences 

The most famous modular form is arguably 

oo 
A(q) =q- 24q 2 + 252q 3 - U72q 4 + 4830<? 5 - . . . = j^ r(n)q n . 

n=l 

Ramanujan discovered a number of congruences involving the coefficients 
r(n), among which is 

r(p) = p +1 (mod 691) for all primes p. 

This is part of a more general phenomenon: if a prime £ > k — 1 divides 
the numerator of the zeta- value £(— k + 1) (equivalently, the numerator of 
the Bernoulli number Bk), then the constant term of the Eisenstein series E^ 
is zero modulo £. This can be interpreted to say that there is a congruence 
mod £ between this Eisenstein series and some cuspidal eigenform of weight 
k. 

As explained in [8] and [9] , Deligne's work on attaching families of £-adic 
Galois representations to Hecke eigenforms of degree one allows us to inter- 
pret Ramanujan's congruence as taking place between traces of Frobenius 
acting on cohomology spaces of local systems. 

A more recent development is the construction (initiated by Laumon 
[181 ] and Taylor [29j] and completed by Weissauer [30(]) of families of four- 
dimensional ^-adic Galois representations attached to Siegel modular eigen- 
forms of degree two. It is then natural to ask about generalizations of 
Ramanujan's congruence to this setting. Building on his study of Eisenstein 
cohomology for arithmetic groups, Harder stated a conjecture [9] regarding 
congruences between classical (degree one) eigenforms and Siegel eigenforms 
of degree two. This statement, which appears below as Conjecture 14.2} was 
verified in a number of cases by Faber and van der Geer [8J, who calculated 
the number of points on the relevant moduli spaces over finite fields and re- 
lated them to the Hecke eigenvalues of Siegel modular forms. (This relation 
was stated as a conjecture in [8(, but it has since been proved by van der 
Geer and Weissauer in most cases. We refer the interested reader to |2J for 
more details.) 

Bergstrom, Faber and van der Geer have extended this point counting ap- 
proach to Siegel modular forms of degree three in [2]. At the same time, they 
formulated another conjectural congruence relating Siegel modular forms of 
degree two and classical eigenforms, this time via critical values of symmet- 
ric square L-functions. We state this below as Conjecture 14.51 It generalizes 
a result of Katsurada and Mizumoto for scalar-valued Siegel modular forms 
(see [13]), and a number of vector- valued cases have been proved by Dum- 
migan in [7, Proposition 4.4]. 

We verify the conjectures by using the data we collected in Section [3] and 
some custom Sage code that computes critical values of L-functions. 
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4.1. Notation. The conjectures appear in different forms in [9|], [8( and [2J. 
We follow the approach of [2] and adapt it to our notation and the quantities 
that we compute. 

We denote the space of cusp forms of weight r with respect to the group 
T^ 1 ' = SL(2,Z) by Sr ■ Suppose / G S r is a Hecke eigenform; we denote 
its Hecke eigenvalue with respect to the operator T(n) by a n = a n (f). The 
spaces Si- are finite dimensional, say of dimension d; according to Maeda's 
conjecture [KJ], the eigenvalues a n are algebraic numbers of degree d. The 
number field that contains the coefficients is denoted Qf. 

Fix a prime p. For an eigenform / e S r , let Qo and ai denote the Satake 
parameters at p and define 

f i pS (f) = a 5 + a 5 a 5 1 for 6 > 1. 

(2) 

Similarly, for a Siegel eigenform F G St, <, let ao, «i, «2 be the Satake 
parameters at p and define 

jy (F) = «q + "o«i + "o«2 + aoOi«2 for 5 > 1. 

4.2. L-functions of modular forms. Let f(q) = Yl a nQ n £ S be an 
eigenform and consider 

oo 

w*)=££=n( i - a *-+* r - i -- j T 1 - 

n=l p 

After introducing the factor at infinity L oa (f, s) = r(s)/(2-7r) s , the completed 
L- function 

A(/, s) = j^pHf, s) = J o f(iy)y a - l dy (4.1) 

has holomorphic continuation to C and satisfies the functional equation 

A(f,s) = (-iy/ 2 A(f,r-s). 

Its critical values occur at 1 < t < r — 1 (of course, the functional equation 
implies that it suffices to consider half of this interval) . 

Manin and Vishik proved that there exist real numbers w_|_(/), w_(/), 
called periods of /, such that the ratio of the critical values of A(/, s) and 
the periods is algebraic. More precisely, define the algebraic critical values 

1 {f t)= f A (f,t)/uj + (f) if t is even 
U ' ; "[A(/,t)/ W _(/) if t is odd. 

Theorem 4.1 (Manin-Vishik[19]). If f £ Sr is an eigenform and t is an 
integer satisfying 1 < t < r — 1, then A(/, t) £ Q/. 

As explained in [9], the denominators of certain classes in the cohomology 
groups of local systems on the moduli space of abelian surfaces should be 
expressed in terms of critical values A(/, t). Harder conjectured that the 
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appearance of certain large primes in these denominators should imply the 
existence of congruences between eigenvalues of forms of degree one and two. 

Conjecture 4.2 (Harder). Let f £ Sr be a Hecke eigenform with coef- 
ficient field Qf and let £ be an ordinary prime in Qf (i.e. such that the 
£-th Hecke eigenvalue of f is not divisible by £). Suppose s € N is such 
that £ s divides the algebraic critical value A(f,t). Then there exists a Hecke 

(2) 

eigenform F 6 SK-, where k = r — t + 2, j = 2t — r — 2, such that 

V (F) = V (/) +p W+i-l) +p*(*-2) (mod £») 
for all prime powers p . 

4.2.1. Computation of the quantities n s. The first step in our numerical 
verification of the conjecture is to compute, as described in Section El the 
Hecke eigenforms in various weights and their corresponding Hecke eigenval- 
ues Ap(-F) and X p 2(F). Once we have those, the second step is to relate the 
Hecke eigenvalues to the values fi p s (F) . We do this by relating the polyno- 
mials that define fi p s(F) to the expressions of X P (F) in terms of the Satake 
parameters ao, a.\, a 2 . 

Example. Let Xi(p 2 ) = Q.{Ti{p 2 ))\ Xi ^ ai . Note X 2 (p 2 ) = p 2fe +^ 6 from the 
definition of the slash operator. Then we have the equations 

X 2 p = X (p 2 ) + ( P + l)X 1 (p 2 ) + (p 2 + l)(p + l)A 2 (p 2 ) 

x p2 = Ac (/) + x 1 ( P 2 ) + x 2 ( P 2 ) 

where the first equation comes from Theorem 12.31 

We compute X p , X p 2 and know X2(p 2 ). This allows us to solve for Xq(p 2 ) 
and X\(p 2 ). Then we observe 

flp — A p 

/i p2 = 2A p2 - A 2 + 2p 2k+ i~ 4 

/V = (3A p2 - 2A 2 + 3(p + l)p 2k+ i- A ) X p . 

4.2.2. Computation of the congruence primes £. We consider the ratios 
A(/,2) : A(/,4) : • • • : A(/,r - 2) and A(/, 1) : A(/, 3) : • • • : A(/,r - 1) 



which by Theorem 14.11 are in Q f . 

Wc compute these ratios of critical values as floating point numbers in 
Sage [25[] . This is done via an implementation of (14. lj) due to Dokchitser [5] . 
We take these floating point numbers and find their minimal polynomial 
using fplll [4], an implementation of the LLL lattice reduction algorithm 
wrapped in Sage. We provide an example to illustrate our process and 
summarize our computations in Table I3~T1 
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Example. Let g G S^ 2 > a two-dimensional space of cusp forms. We will 
calculate the ratio of critical values 

A(g, 1) : A(5, 3) : A(g, 5) : A(g, 7) : A( 5 , 9) : A(g, 11) : A( 5 , 13) : A(g, 15). 

As we are interested only in the ratio, we compute 

A(g,l) . A(ff,3) . A(g,5) . A(g,7) . A(g,9) . A(g,ll) . A(g,13) . A(g,15) _ 
A(g,l) • A(g,l) " A(g,l) ' A(g,l) ' A(g,l) • A(g,l) ' A(g,l) ' A(g,l) ~ 

1 : 0.045375 • • • : 0.002369 • • • : 0.000143 • • • : 0.000010 . . . 

8.65221 • • • x 10~ 7 : 8.50052 • • • x 10~ 8 : 9.23745 x 10~ 9 

using a Sage implementation of A(g, s). Then we find the minimal polyno- 
mial of each ratio; e.g., -AU has minimal polynomial 

1254224510x 2 - 471820065x + 18826702. 

4.2.3. Checking ordinarity of the primes £. We use a simple algorithm that 
is very fast but uses large amounts of storage. Let d be the dimension of Sy. , 
and suppose we want to check that £ is ordinary for all Hecke eigenforms in 
Sr ■ We proceed as follows: (a) compute the Victor Miller basis for Sf 
to a precision of about d£ coefficients; (b) compute the matrix of the Hecke 
operator Ti acting on this basis; (c) reduce the matrix modulo £ and check 
whether it is invertible. 

This allowed us to verify that most primes £ appearing in Table 14.11 are 
ordinary. The current understanding of the distribution of non-ordinary 
primes is rather limited, but numerical evidence seems to indicate that they 
are very rare in level one, so it would be surprising to find a prime £ > r 
that divides an algebraic critical value and is non-ordinary 

4.2.4. Verification of the congruences. We observe that there are two cases: 
when (5 = 1 and when 6 > 1. The difference is that in the first case the 
congruence reduces to a congruence on the Hecke eigenvalue A p while in the 
second case we require both X p and X p 2 . The effect of this difference is that 
we can verify many more congruences when 5 = 1 than when 6 > 1; this is 
due to the number of coefficients needed to compute \ p as compared to the 
number needed to compute X p and X p 2. 

The way the actual verification works is essentially the same starting from 
the point where n p s(F) and /J,„s(f) have been computed. Each side of the 
congruence mod £ s is an algebraic number in <Qp and Q/ respectively. In 
the cases we have considered, the exponent s appearing in the conjecture 
was always 1. We compute the minimal polynomial m(x) of the coefficient 
a p of / and the minimal polynomial M(x) of the Hecke eigenvalue X p of F. 
Then we look at the roots of m and M in F^. The conjecture holds if for 
some choice of root of m and some choice of root of M the congruence holds. 

The following statement summarizes our results on Conjecture 14.21 
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Theorem 4.3. Let r < 60 be a multiple of 4. If f € Sr is a Hecke 
eigenform with coefficient field Qf and £ is an ordinary prime in Qj that 
divides the algebraic critical value A(/, r/2 + 2), then there exists a Hecke 

(2) 

eigenform F £ S A 2 such that 



for 



/V (F) ee „ pS (/) + ^/ 2+ D + /(''/2-2) (mod 



p 5 e {2, 3, 4, 5, 7, 8, 9, 11, 13, 17, 19, 23, 25, 27, 29, 31, 125}. 



Proof. For weights r < 28, we have verified that there are no ordinary primes 
I dividing A(/, r/2 + 2), so the statement is vacuously true. 
For weights 32 < r < 60 and 

p s G {2, 3, 4, 5, 7, 9, 11, 13, 17, 19, 23, 25, 29, 31} 

we have verified the congruence for all large primes £ dividing the algebraic 
critical value. The results are listed in Table I4~T1 

The remaining cases p € {8, 27, 125} follow from the rest by Proposi- 
tion KB □ 



32 



36 



40 



44 



48 



52 



56 



60 



18 



20 



22 



24 



2<i 



28 



30 



32 



large £ \ Norm(A(/,t)) 



211 



269741 



509 

1447 



205157 



168943 



173 

929 

4261 

* 434167 



173 

1721 
38053 

1547453 



* 325187 

* 32210303 

* 427092920047 



(k,j) 



(16,2) 



(18,2) 



(20,2) 



(22,2) 



(24,2) 



(26,2) 



(28,2) 



(30,2) 



dimS' 



W 

k,j 



10 



11 



Table 4.1. A summary of the cases verified numerically for 
the proof of Theorem 14.31 (The primes £ marked with a * 
have not been checked to be ordinary.) 
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4.3. Symmetric square L-functions of modular forms. It is possible 
to associate higher-degree L- functions to modular forms, by using various 
tensorial constructions. We describe the L-function attached to the sym- 
metric square of a modular form. 

Fix a Hecke eigenform / G Sr and let a p , (3 P be the roots of the polyno- 
mial X 2 — a p X + p r ~ l . The associated symmetric square L-function is 

L(Sym 2 /, s) = H ((1 - a 2 |T s )(l - /? p V s )(l - Op&p-)) -1 . 

v 

We take as factor at infinity 

L (gynrV,)- r(*)r(( a + 2-r)/2) 

and set 

A(Sym 2 /, s) = L^Sym 2 /, s)L(Sym 2 f, s). 

Then A(Sym /, s) has holomorphic continuation to C and satisfies the func- 
tional equation 

A(Sym 2 /, a) = A(Sym 2 /, 2r - 1 - s). 

We define the algebraic critical values 

A(Sym 2 /, t) = ff^4 for t = r, r + 2, . . . , It - 2, (4.2) 

where (•, •} denotes the Petersson inner product. (It is possible to ex- 
press the algebraic critical values as quotients of the completed L-function 
A(Sym /, t), which would be closer to the treatment of the usual L-function 
as given in the previous section. We prefer to take a quotient of L(Sym 2 /, t) 
instead, as this is the definition used in much of the existing work on sym- 
metric square L- values.) 

Theorem 4.4 (Zagier [3l|). If f € Sr is an eigenform and t is even such 
that r < t <2r — 2, then A(Sym /, t) is an algebraic number. 

Moreover, and this is useful for our computations, it can be shown that 
A(Sym 2 /,t)GQ / , see 0. 

Conjecture 4.5 (Bergstrom-Faber-van der Geer-Harder). Let f € Si- be 
a Hecke eigenform with coefficient field Qj and let £ be a large prime in Qj . 
Suppose s € N is such that £ s divides the algebraic critical value A(Sym /, t). 

(2) 

Then there exists a Hecke eigenform F S S^, ■, where k = t — r + 2, j = 
2r — t — 2, such that 

/VtL^^fX/^+l) (modf) 

for all prime powers p s . 
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The case j = concerns scalar-valued Siegel modular forms. The first 
examples of such congruences were found by Kurokawa [la], who conjectured 
that they should be governed by certain primes dividing the numerators 
of algebraic critical values. Kurokawa's conjecture was recently proved by 
Katsurada and Mizumoto, who even extended these results to the case of 
scalar-valued Siegel modular forms of arbitrary degree (see [la . Theorem 
3.1]). 

In the vector- valued setting, the congruence in Conjecture 14.51 was proved 
for the six rational eigenforms of degree one (weights 12, 16, 18, 20, 22, 
26) by Dummigan in [7, Proposition 4.4]. (Dummigan has indicated that it 
should be possible to extend his Proposition 4.4 to higher weights, using a 
pullback formula as in Katsurada and Mizumoto [l3j.) 

Remark 4.6. The conjecture does not specify what is meant by a large prime 
I. Dummigan's result uses £ > 2r. In the cases we have verified (see Theo- 
rem H21 for details), it was sufficient to take £ > 2. 

Our numerical verification of Conjecture 14.51 follows the approach of the 
last section. We highlight only the essential differences. 

4.3.1. Computation of the congruence primes £. We find the appropriate 
primes £ by computing the algebraic critical values directly from (14. 2p . The 
squared-norm (/, /) of / can be obtained from the identity 

(/,/) = ^fe^T^(Sym 2 /,r), 

so all we require is high-precision evaluation of the symmetric square L- 
function at various points. For this we use Dokchitser's L-function calcula- 
tor g as wrapped in Sage, as well as some Sage code made available to us 
by Martin Raum. 

Having obtained a sufficiently precise floating point approximation to the 
algebraic number A(Sym f,t), we then find its minimal polynomial. The 
congruence primes £ are the primes larger than 2 occurring in the factoriza- 
tion of the numerator of the norm of A(Sym /, t). 

The critical values we obtain in this way agree with the ones computed by 
Dummigan in the case of rational eigenformqj, see Table 1 in [6y . We were 
also able to verify the case r = 24 by comparing our result with the trace of 



A(Sym /, 46) as obtained (by theoretical means) by Lanphier in 17] 



The following statement summarizes our results on Conjecture 14.5 



Dummigan confirmed that a few of the factorizations from Table 1 in [6| are incorrect. 
Here are the values in question, with their corrected factorizations: 

fc = 16,r = 3: 2 20 /3 7 • 5 3 • 7 ■ 11 ■ 13 2 ■ 17 

k = 16, r = 11 : 2 24 • 839/3 12 • 5 8 • 7 4 • ll 2 • 13 2 • 17 • 19 • 23 

k = 20, r = 11 : 2 27 • 304477/3 19 ■ 5 
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Theorem 4.7. Let r < 32. If f G Sr is a Hecke eigenform with coefficient 
field Qf and I > 2 is a prime in Qf that divides the algebraic critical value 

A(Sym /, 2r — 4), then there exists a Hecke eigenform F G S^_ 2 2 such that 



for 



V (F)=/v(/)(/ (r - 4) + l) (mod£) 



p 5 G {2, 3, 4, 5, 7, 8, 9, 11, 13, 17, 19, 23, 25, 27, 29, 31, 125}. 



Proof. For weight r = 12, we computed the numerator of the rational num- 
ber A(Sym A, 20) and found it to be — 2 23 , so there are no large primes 
dividing this algebraic critical value and the statement is vacuously true. 
For weights 16 < r < 32 and 

p s G {2, 3, 4, 5, 7, 9, 11, 13, 17, 19, 23, 25, 29, 31} 

we have verified the congruence for all large primes £ dividing the algebraic 
critical value. The results are listed in Table 1431 

The remaining cases p G {8, 27, 125} follow from the rest by Proposi- 
tion EEEJ □ 



r 


t 


odd£ | Norm(A(Sym 2 /,t)) 


(k,j) 


dimSg 


16 


28 


373 


(14,2) 


1 


18 


32 


541 
2879 


(16,2) 


2 


20 


36 


439367 


(18,2) 


2 


22 


40 


281 
286397 


(20, 2) 


3 


24 


44 


2795437 
256021114049 


(22,2) 


5 


26 


48 


4598642018203 


(24, 2) 


5 


28 


52 


4017569791 
65593901428085768723 


(26,2) 


8 


30 


56 


937481 
4302719815755987715030485446839 


(28,2) 


10 


32 


60 


350747 

45130901953 

432796809552670722149 


(30, 2) 


11 



Table 4.2. A summary of the cases verified numerically for 
the proof of Theorem 14.71 
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4.4. Reduction of cubes to primes and squares of primes. We de- 
scribe some elementary considerations that allow reducing the case 5 = 3 of 
both conjectures to the cases 5 = 1 and 5 = 2. For ease of notation in this 
section, we will write 

95 = «o + «o a i 
G$ = a + a a 1 + a a 2 + a a 1 a 2 , 

where in the first line «o an d (X\ are the Satake parameters of / G <Sf , while 

(2) 

in the second line ao, «i and a 2 are the Satake parameters of F € S k ■. 

In the degree one setting, we have the relation a^a\ = p r ~ l , which allows 
us to express g 2 and 53 in terms of g\ : 

g 2 = g\ - 2f-\ g 3 = 9l (g 2 - Sp^ 1 ). (4.3) 

In the degree two setting, we have the relation o^a\(X2 = p 2k+ i~ 3 , which 
allows us to express G3 in terms of Gi and G 2 : 

G 3 = \d (-G\ + 3G 2 + 6p 2fc+ ^ 3 ) . (4.4) 



Proposition 4.8. In Conjecture \4-%\ and Conjecture \4-5j the congruences 
for the case 5 = 3 follow from the congruences for the cases 5 = 1 and 5 = 2. 

Proof. 

(a) Define h$ = p^+i -1 ) + p°( k ~ 2 > for all 5 > 1. Then the congruence 
in Conjecture 14.21 can be written 

(CO: Gs = gs + h 5 (mod f). 

It is easily seen that 

h 2 = h\- 2p 2k+j - 3 , h 3 = h x (h\ - 3p 2k+j - 3 ). (4.5) 

(Observe the similarities between these equations and ()4.3p ,) 

We assume that the congruences (Ci) and (C2) hold. Using Equa- 
tions (14. 4p . (I4.3J) and (14.51) (in this order), we compute 

G 3 = \{gi + h x ) (-{ 91 + hrf + 3( 52 + h 2 ) + Qp 2k+ ^\ 
= o(5l + hi) {{-9 2 i + 352) - 2g 1 h 1 + {-h 2 + 3/i 2 + Qp 2 



2 

= {gi + h 1 ){g 2 -3p r - l -g 1 h 1 + h 2 ) 

= 9l- 3p P-1 <7i + h\- 3p p - 1 /n 

= 93 + h 3 , 

after noting that, under the conditions of Conjecture 14.21 the weight 
parameters r, k and j are related by r — 1 = 2k + j — 3. 
(b) The calculation is similar to the previous part. We let h$ = p 5 ( fc-2 ) + l 
for all 5 > 1. The congruence in Conjecture 14.51 takes the form 

(C' s ): G s = g 5 h s (mod^ 2 ). 
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We easily see that 

h 2 = h\-2p k ~ 2 , h 3 = h 1 ( y hi-3p k - 2 ). 
Assuming that congruences (C[) and {C' 2 ) hold, we obtain 

G 3 = ^i (-^ + 3( 5 ? - 2/- 1 )(/ i ? - 2p k ~ 2 ) + 6p 2fc+ ^ 3 

= gihx {glh\ - 3p k ~Vi ~ W- X h\ + 9p 2k+ ^ 3 

= 9ihi(£-3p r - 1 )(h 2 -3 P k - 2 ) 

= 53/13, 

where we used the relation r = k + j, valid under the conditions of 
Conjecture 14.51 

□ 
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